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Abstract 

The development of flux-vector splitting through the 1970s and 1980s is re- 
viewed. Attention is given to the diffusive nature of flux- vector splitting, 
which makes it an undesirable technique for approximating the inviscid fluxes 
in a Navier-Stokes solver. Several proposed improvements, including a brand- 
new one, are discussed and illustrated by a simple, yet revealing, numerical 
test case. Finally, an outlook for flux- vector splitting in the 1990s is pre- 
sented. 


1 Introduction 

Flux-vector splitting is a technique for achieving upwind bias in numerical 
flux functions for the Euler equations. It has been very popular during the 
Euler era (the 1980s), but has lost much of its appeal with the rise of Navier- 
Stokes solvers, since it causes undesirably high errors in boundary layers. 
Ironically, it shares this deficiency with other flux functions that are widely 
used in Navier-Stokes codes [1]. 

Flux- vector splitting is a natural consequence of regarding a fluid as an 
ensemble of particles. Measured along any coordinate, some particles will 
move forward, other ones backward; this automatially splits the fluxes of 
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mass, momentum and energy into forward and backward fluxes, i.e. 


F = F + + F" 


( 1 ) 


On a computational grid these split fluxes immediately allow upwind differ- 
encing: 
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( 2 ) 


A full update step with the first-order upwind flux-split scheme for the Euler 
equations can be regarded as an approximate way to integrate the collisionless 
Boltzmann equation 


d±(q) + , d(j)(g) = Q 


( 3 ) 


dt Sx 

for the distribution function <f>(q ) of the particle velocity q. For the sake 
of numerical simplicity the velocity distribution in practice usually is not 
chosen to be Maxwellian (see, however, [2]). Because of this, and the finite 
length of the time-step, the velocity distribution in any cell at the end of a 
time-step is not of the same functional form as at the start of the time-step, 
when it was assumed to be an “equilibrium” distribution. Replacing the 
final, non-equilibrium distribution by an equilibrium distribution with the 
same integrated mass, momentum and energy in each cell instantaneously 
simulates the effect of the particle collisions. Euler schemes based on flux- 
vector splitting may therefore be termed Boltzmann- type schemes [3]. 


2 A gallery of flux- vector splittings 

2.1 The Beam Scheme, alias Steger- Warming split- 
ting 

The first Boltzmann-type scheme in use for integrating the Euler equations 
was developed for astrophysical calculations around 1970 by Kevin Prender- 
gast [4]; it is based on a velocity distribution consisting of delta- functions, 
whence the name “Beam Scheme”: 

M = pA<1 -{ u - o-)] + ( p - 2 p s )8(q -u) + p s 8[q - (u + cr)]. (4) 


204 



Here p is the density, p s is the density of the particles in each side beam, u is 
the mass-averaged flow velocity, and a is a dispersion velocity not less than 
the speed of sound. In the special case that a equals the speed of sound, 
and the specific internal energy that is not in the form of ^-translations (in 
the Beam Scheme a passively convected quantity) is equidistributed over the 
side beams, the flux splitting becomes identical to the well-known splitting 
of Steger and Warming [5], developed in the late 1970s without regard to a 
velocity distribution. This equivalence was discovered by Van Albada [6]. 

The Beam-Scheme interpretation is advantageous, for instance, when for- 
mulating the flux splitting for gases with an equation of state more compli- 
cated that the ideal-gas law [7]. 

2.2 Van Leer’s splitting 

The most popular splitting subsequently developed, Van Leer’s [8], again is 
not derived from or associated with any a priori chosen velocity distribution, 
but from certain mathematical constraints. The main design features of Van 
Leer’s splitting are: 

® the split fluxes are continuously differentiable (Beam/ Steger- Warming 
fluxes are non-differentiable in sonic and stagnation points, which is 
inappropriate when approximating a first-order system of conservation 
laws); 

m for subsonic flow the Jacobians of the split fluxes have a zero eigenvalue, 
which accounts for crisp numerical profiles of stationary shocks. 

If the second constraint is relaxed, a one-parameter family of continuously 
differentiable split fluxes can be generated [9]; these are the simplest possible 
in the sense that they are at most quartic in the Mach number, just as the 
Van Leer fluxes. The differences among members of this family arise only in 
the energy-flux splitting. 

Extensions of these fluxes to real gases are also included in [9], but have 
been given by several other authors [7] [10]. 
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2.3 HanePs energy-flux splitting 

The above family of energy-flux splittings includes the one originally pro- 
posed by Hanel [11]: 

Eli Ei± T I 

^energy 1 mass 1 1 1 

where H is the specific total enthalpy. Advantages of this formula are: 

• it is as simple as can be; 

• it admits steady Euler solutions with constant total enthalpy through- 
out the flow. 

Hanel claims that this flux splitting, when used in Navier-Stokes calcula- 
tions, gives more accurate total-enthalpy values in the boundary layer. This 
may have been observed for the lower flow speeds; in the hypersonic flow 
regime the improvement is insignificant (see further in Section 5). 

3 Flux-vector splitting as an approximate 
Riemann solver 

Flux-vector splitting can be used beyond the first-order upwind scheme; in 
a higher-order MUSCL-type code [12] it simply takes the same place as any 
“approximate Riemann solver,” such as Roe’s [13] or Osher’s [14]. It is one 
more way to merge the two state vectors on the left and right side of a cell 
interface into one flux vector, namely, by 

F(U £ ,U,i) = F+(U t ) + F-(U*). (6) 

4 Dissatisfaction with flux- vector splitting 

Dissatisfaction with flux-vector splitting in Navier-Stokes codes stems from 
the dissipative properties of flux-vector splitting: contact discontinuities and 
slip surfaces are diffused, even when these are stationary and aligned with 
the grid [8] [3]. This is in contrast to flux-difference splittings such as Roe’s 
and Osher’s. A serious consequence is that the attached boundary layers to 
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be resolved by Navier-Stokes codes get artificially broadened, and adiabatic- 
wall temperatures become inaccurate; numerical solutions obtained with flux- 
difference splitting are vastly superior [1]. 

This can be easily understood from, for instance, the net transverse- 
momentum flux across the boundary layer, as computed with flux-vector 
splitting ( B denotes bottom cell, T top cell): 


pvu : = F^ b u b + F massT u T 

= (FL..B + 


(T) 


( 8 ) 


When the net mass flux 


rpnet 
* mass 


-^mass B T ^mass X 


(9) 


vanishes, the second term in the momentum flux still causes numerical dif- 
fusion. 


5 Improvements to flux formulas 

Hanel [15] has suggested to replace the above formula for the transverse- 
momentum flux by one borrowed from flux- difference splitting: 


pvu := 

-^mass^upwind ) 


(10) 

^upwind 

U B 

if J?net 
mass 

>o, 

(11) 

^upwind — 

UT 

jf 

mass 

< 0 . 

(12) 


This mixture of flux-vector splitting and flux- difference splitting prevents 
the numerical broadening of the boundary layer, but does not improve the 
accuracy of the wall temperature. It further introduces pressure irregularities 
across the boundary layer. 

Nevertheless, the partial success of Hanel’s modified transverse- momentum 
flux suggests that a further improvement can be obtained by introducing a 
similar formula for the energy flux, i.e.: 

pvH := F^f s H upwhvh (13) 
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with 


-^upwind 

H b 

if p net 

mass 

>0, 

(14) 

^upwind 

H t 

if p net 
mass 

< 0 . 

(15) 


A numerical test confirms the improvement, as shown in Figures 1 and 2. 
Plotted are the temperature and pressure distributions from the solution to 
the conical- flow problem of [1], as computed with the Roe fluxes, the Van 
Leer fluxes, the two Hanel modifications of Equations 5 and 10, and the extra 
modification 13 suggested in this section. The wall temperature computed 
with Equation 13 is significantly closer to the correct value of 11.7 than for 
the other Van Leer-type splittings; unfortunately, the pressure irregularity 
remains. Convergence histories for the four calculations are shown in Figure 
3 and indicate that the modified schemes are as robust as the original Van 
Leer scheme. 

It is clear that we have hardly begun to investigate mixtures of flux- 
vector and flux-difference splitting; the success of the formula introduced in 
this paper suggests that there is something to gain here. The advantage of a 
simple blended formula over full flux-difference splitting is that it will allow 
a complete linearization for use in implicit schemes; for Roe’s flux formula, 
for instance, only an incomplete linearization is acceptable in practice (see 
the horrendous complete formula in [16]). 

6 Can pure flux-vector splitting be saved? 

However promising the above improvement may be, an intriguing question 
still remains: 

• Is it possible at all to construct a flux-vector splitting that does not 
diffuse a grid-aligned boundary layer? 

This question can be rephrased as: 

• Is it possible to split the Euler fluxes such that both F+ ass and 
vanish with the flow speed? 
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a Roe 
s Van Leer 
A Hanel 89 
„ VL/H 90 


Figure 1: Numerical solution of the Navier-Stokes equations for self-similar 
hypersonic flow (Moo = 7.95) over a 10° cone; shown are the non-dimensional 
temperature distributions computed with four different flux functions. 
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Figure 2: Same solutions as in Figure 1; pressure distributions. 
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Conical flow, first-order scheme 



Iterations 


Figure 3: Residual histories of the solutions in Figures 1 and 2. 
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If indeed this were possible, the form of the flux splitting for small v 
would follow immediately from symmetry considerations: 


F ± = 

mass 

\pv + 0(v 2 ), 

(16) 

Z7»± 

mom|| 

± apcv + 0(v 2 ), 

(17) 

Tp± 

1 moml 

puv + 0(v 2 ), 

(18) 

F ± = 

energy 

i pHv + 0(v 2 ), 

(19) 


where c is the sound speed and a. is a positive-valued free parameter. For 
vanishing a this splitting leads to central differencing, which will be neutrally 
stable for v — 0 and unstable for v ^ 0, if foi'ward time-differencing is used. 
Positive values of a will introduce some dissipation, but it is clear that at 
least a multi-stage time-stepping algorithm must be used to render stability. 

A greater problem is that the above flux values for small v must smoothly 
join the branches for larger v, that is, the standard Van Leer-Hanel fluxes. So 
far, the numerical experience is that the overall flux splitting can not be sta- 
bilized. This result, however, is not conclusive; the possibility of constructing 
a continuously differentiable flux splitting that reduces to the above formu- 
las for slow flow, and is stable for all flow speeds when used in a suitable 
time-marching scheme, still remains open. 

7 Multi-dimensional flux-vector splitting 

The interpretation of flux-vector splitting as a consequence of positive and 
negative particle speeds included in a distribution function makes it possible 
to extend the concept to the multi-dimensional Euler equations in a grid- 
independent way. For instance, when extending the Beam Scheme into two 
dimensions, we may consider introducing just one “beam” moving with the 
average flow velocity, and a circular “front” moving, relative to the average 
flow, with the dispersion velocity: 

^(9®) Py) = Pb eam^ u ) 2 + (?» - vY^j+PirontS (y/(qx ~ u) 2 + (q y ~ v) 2 - a 

( 20 ) 


211 



The standard two-dimensional extension of the distribution function calls 
for five beams; the two added beams represent the velocity dispersion in 
the ^-direction. This introduces a directional bias in numerical solutions; in 
particular, shock waves oblique to the grid may be excessively smeared. The 
above distribution function should eliminate this effect; the penalty for its 
use is that the expressions for the fluxes out of a cell are more complicated. 

Extension of the Van Leer fluxes in a grid-independent way is more prob- 
lematic, as these do not derive from a given distribution function. It is 
possible, a posteriori , to find some velocity distribution that explains the 
one-dimensional Van Leer fluxes; this could then be extended in an omni- 
directional way. The resulting fluxes, however, would be somewhat arbitrary, 
and very complicated. It would be preferable to derive the multi-dimensional 
version from purely algebraic considerations; this approach remains to be in- 
vestigated. 

8 An outlook for flux- vector splitting 

Flux-vector splitting is still alive, even for Navier-Stokes applications. It may 
be mixed with flux-difference splitting, for the sake of preventing numerical 
diffusion of boundary layers, and there may still be a way to avoid this 
diffusion in a pure flux-vector splitting. More work needs to be done; the 
reward could be a simple and robust, easily linearized inviscid-flux formula, 
with an accuracy rivalling that of Roe’s flux formula. 
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